Competition between Antiferromagnetism and Superconductivity in High T c Cuprates. 
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Using variational cluster perturbation theory we study the competition between d-wave super- 
conductivity (dSC) and antiferromagnetism (AF) in the the t-t'-t"-U Hubbard model. Large scale 
computer calculations reproduce the overall ground state phase diagram of the high-temperature 
superconductors as well as the one-particle excitation spectra for both hole- and electron-doping. 
We identify clear signatures of the Mott gap as well as of AF and of dSC that should be observable 
in photoemission experiments. 



Ever since the first paper of P.W. Anderson on high- 
temperature superconductivity [JJ, most theorists have 
been working with the premise that the physics of these 
most intriguing materials should be contained in the 
two-dimensional one-band Hubbard model. That model 
of strongly interacting electrons on a two-dimensional 
square lattice contains a kinetic energy term that rep- 
resents the band structure, and a potential energy term 
that accounts only for on-site (perfectly screened) inter- 
action. This simplicity may be too naive, in particular 
in light of the striking contrast between the phase dia- 
gram of the electron-doped and that of the hole-doped 
cuprates. The former has an antiferromagnetic (AF) 
phase that extends over a wide range of doping while the 
d-wave superconducting (dSC) phase seems restricted to 
a small range. The reverse is true for the hole-doped 
case. In addition, some hole-doped cuprates seem to ex- 
hibit stripes @ or other recently discovered phases Q 
near the boundary between antiferromagnetism and d- 
wave superconductivity. We cannot trust that copper- 
oxygen planes described by the Hubbard model contain 
the whole physics until we have shown that both the 
hole- and electron-doped cuprates can be explained by 
this model. Yet previous theoretical studies that found 
agreement with experiment concentrated on the hole- 
doped cuprates. In addition, many authors argue that 
the one-band Hubbard model description of the cuprates 
is inadequate to obtain dSC |aJg|or that additional phys- 
ical effects must be invoked pjlal- 

In this letter we show, in agreement with the basic 
premise of Anderson, that the two-dimensional Hubbard 
model accounts for the observed doping dependence of 
the AF and dSC phases as well as for the single-particle 
excitations for both electron and hole doping. The strik- 
ing asymmetry between electron and hole doping can be 
accounted for by the fact that the known band structure 
is not particle-hole symmetric and that interactions at 
intermediate-coupling cannot completely wipe out band 
effects. We also identify distinct spectral features of 
Mott, AF and dSC physics. 

Our results were obtained thanks to a new method- 
ology, namely Variational Cluster-Perturbation Theory 
(V-CPT) 0,' and because of remarkable increases in com- 



puter power made possible by cluster architectures. We 
study the Hubbard model with interaction U and with 
hopping parameters that are taken from band structure 
calculations namely a diagonal hopping t = — 0.3i 
and third-neighbor hopping t" — 0.2t. 

Variational cluster perturbation theory. V-CPT is an 
extension of Cluster Perturbation Theory (CPT) 
that is based on the self-energy-functional approach 
(SFA) [3. This approach uses the rigorous variational 
principle <5ftt[S]/5£ = for the thermodynamic grand- 
potential fit written as a functional of the self-energy £ 



Qt[E] = J F[S]+Trln(-(Go 1 -E)- 1 ). 



(1) 



The index t denotes the explicit dependence of fit on 
the matrix tjj of hopping terms or, more generally, on all 
one-body operators. That dependence comes through the 
Green function Go of the one-body part of the Hamil- 
tonian. In the above expression, F[E] is a universal 
functional of the self-energy obtained from the Legen- 
dre transform of the Luttinger-Ward functional. The 
physical Green function is G = —5F/8Yi and the sta- 
tionary condition for f2t[E] give Dyson's equation. Al- 
though F is universal, its exact form is unknown. But 
all Hamiltonians with the same interacting part share 
the same functional form of F[T]. Hence F[S\ may be 
derived from the exact solution of a simpler Hamiltonian 
H' whose choice of one-body terms makes it exactly solv- 
able. One then looks for stationary solutions within the 
subspace defined by that simpler solvable problem. The 
Dynamical Mean-Field theory (DMFT) [D| is a special 
case [2( where the simpler problem (H') is that of a sin- 
gle site (the "impurity" ) connected to an infinite bath 
of non-interacting sites. When the simpler problem con- 
sists of a single cluster connected to a bath, one recov- 
ers Cellular-DMFT Il4l . The Dynamic Cluster Approx- 
imation (DC A) 0, llfil does not fit in that general 
variational scheme although it is self-consistent. In 
V-CPT, one uses for H' a Hamiltonian formed of clus- 
ters that are disconnected by removing hopping terms 
between identical clusters that tile the infinite lattice. In 
the cluster, the Hamiltonian contains the original Hub- 
bard interaction, hopping, and Weiss fields that are go- 
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ing to be determined by minimizing the grand potential 
Sit- In V-CPT the cluster shape and size can be varied 
as a check of the validity of the results in the thermo- 
dynamic limit. V-CPT is thcrmodynamically consistent 
and causal. It goes beyond ordinary mean-field theory, 
since it predicts the absence of long-range AF order in 
one dimension 

In practice, one proceeds as follows. The universality 
of F allows us to express the functional fit in terms of the 
corresponding functional for the cluster Hamiltonian H', 
i.e., flf in Eq.(5) of Ref. 0. To include the possibility 
of broken symmetry, H' contains the terms 

H M = M£(-l) (na T -n a ;) (2) 

a 

H D = A «fc ( c aTCfci + H.c.) (3) 

a , 6 

where A a & = D if sites a and b are nearest-neighbors 
(NN) along the x axis and A a & = — D if the sites are NN 
along the y axis. Such a Weiss field cannot be obtained 
from a Hartree-Fock factorization of the Hubbard model. 
The Weiss fields M for AF and D for dSC phases allow 
for the physics of long-range order while the dynamical 
short-range interactions are taken into account in the ex- 
act diagonalization of the cluster. One then restricts the 
space of self-energies to the exact self-energies of all pos- 
sible cluster Hamiltonians: E = S(t'). In that case the 
set of variables t' includes the variational parameters M 
and D. Rearranging Eq.(5) of Ref. 12] in terms of the 
cluster grand-potential, Q! = f2[[E] and Green function 
G' -1 = Gq -1 - E, one finds 

fMt')=W" ^$>det(l + (Go 1 - GfT 1 ^'), (4) 
Jc K 

which is the starting point of numerical calculations. The 
functional trace has now become an integral over the di- 
agonal variables (frequency and superlattice wave vec- 
tors) of the logarithm of a determinant over intra-cluster 
indices. The frequency integral is carried along the imag- 
inary axis, since it can be shown that the integrand 
decreases asymptotically as l/(iuj) 2 . The minimum of 
fi(M, D) is found using the conjugate-gradient algorithm. 
Computing time varies a lot, from one cpu-hour in the 
case of a 4-site cluster, to hundreds of cpu-hours for a 
10-site cluster. 

Order parameters. The treatment of superconductiv- 
ity in V-CPT is best done in the Nambu representation. 
Translation invariance in the superlattice implies that the 
Green function G _1 = Go -1 — E depends, after Fourier 
transformation, on two wave vectors k and k' that can 
only differ by a reciprocal superlattice vector. The elec- 
tron concentration (n), the AF order parameter (M ) and 
the dSC order parameter (D a ) are respectively expressed 



as 

f dk f duj sr^ „ . 
n = y c -X>(k,k,.) (5) 

where Q = (n, n) is the antiferromagnetic wave vector, 
Q a the normal and T the anomalous V-CPT Green func- 
tions, written now as functions of wavevector. The dSC 
form factor is g(k) = cos(k x ) —cos(k y ). The integral over 
frequencies is carried along the upper half of a clockwise 
contour that encloses the poles of the Green function up 
to uo = 0. 

Results. We present V-CPT calculations on L = 6- 
site (2 x 3) , 8-site (2 x 4) clusters and the 10-site cluster 
of Ref. [17| , using M and D as simultaneous variational 
parameters and U and /i (the chemical potential) as con- 
trol parameters. The number of electrons within the 
cluster (N c ) is a conserved quantity only when the dSC 
Weiss field D vanishes. When D ^ 0, n c (n c — N c /L) 
can take any real value controlled by //. The calculated 
electron concentration n (Eq. [SJ is not quantized, even 
when D — 0, except when /x lies within the Mott gap 
(n = n c = 1). Even though results are displayed as func- 
tions of n, one must keep in mind that it is the chemical 
potential fj, that is the true control parameter. 

The Weiss fields M and I? as a function of n are ob- 
tained from the stationary point (in this case a mini- 
mum) of the grand potential. The Weiss fields gener- 
ally decrease with increasing cluster size, as it should for 
spontaneous symmetry breaking. Instead of the Weiss 
fields, we show the order parameters as a function of the 
electron concentration n calculated from Eq. [SI for dif- 
ferent system sizes in Fig. Q and for different interaction 
strengths in Fig. [21 

The bottom panel showing Mq in Fig. ^ clearly repro- 
duces the striking contrast between electron- (n > 1) and 
hole- (n < 1) doped systems. Antiferromagnetism can 
extend up to about 15% doping on the electron-doped 
side while it persists up to only about 6% on the hole- 
doped side. The dependence on system size is not mono- 
tonic because even at fixed cluster size there can be a 
dependence on the shape of the cluster. Nevertheless, 
there is a clear convergence with system size, especially 
on the electron-doped side. The dSC order parameter 
is shown on the top panel. For different system sizes, 
the vertical lines indicate the filling where the dSC phase 
appears by itself, without AF order parameter. In the 
electron-doped case, the pure dSC phase is on the right. 
It survives for a narrow range only, namely n = 1.13 
to n = 1.15 for L = 8 for example. In the hole-doped 
case, the dSC phase on the left of the vertical lines is 
present for a range of fillings that is at least twice as 
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FIG. 1: AF (bottom) and dSC (top) order parameters for 
U = 8t as a function of the electron density (n) for 2x3, 
2x4 and 10-site clusters. Vertical lines indicate the first 
doping where only dSC order is non- vanishing. 
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FIG. 2: AF and dSC order parameters as a function of the 
electron density on a 2 x 4 cluster for U = 6t, U = 8t and 
U = 12*. 



large, n = 0.87 to n = 0.93 for L = 10 for example. Un- 
fortunately, finite cluster effects do not allow us to obtain 
reliable results for larger dopings for either the hole or 
electron cases. The dSC order parameter is not present 
on bonds between clusters, making the results more sen- 
sitive to size effects than in the case of the AF phase. 
Fig- also shows, in the electron-doped AF+dSC 
phase where AF and dSC order parameters are both non- 
vanishing 0,0, E3 We verified that, as expected from 
symmetry, the 7r -trip let order parameter is non- vanishing 
in that phase |2ll,l22| . which is separated from a pure dSC 
phase by a quantum critical point around 13% doping, 



near the value suggested by experiment |23j, |2J] . The sit- 
uation is less clear on the hole-doped side where the L = 6 
cluster has a very small doping range for the AF+dSC 
phase, the L = 8 cluster a large one, while the L = 10 
cluster shows none. This suggests that the way in which 
the AF and dSC phases approach each other on the hole- 
doped side cannot be accurately described by the small 
variational space that we use. Additional order parame- 
ters, such as stripe Q or checkerboard orders Q observed 
in certain cuprates may be necessary to get the full pic- 
ture. No SO(5) symmetric point ;25] appears in our cal- 
culation in a size-independent way. On the other hand, 
our results for Dq in Fig. ^ show unambiguously that 
the pure dSC phase appears over a much broader range 
of dopings for hole- than for electron-doped cuprates, as 
observed experimentally. 

It is also instructive to know how the ground-state or- 
der parameters vary with interaction strength U , espe- 
cially because several normal-state calculations for the 
pseudogap [i^, I22} show that the interaction strength for 
electron-doped cuprates near optimal doping should be in 
the weak to intermediate coupling range (U ~ 6t), with 
U increasing as n decreases. A look at Fig. |2 for D and 
M shows that the range of dopings where only D is 
non-vanishing is larger on the hole than on the electron- 
doped side for all values of U. That range increases with 
U in all cases so that a drop in U as n increases rein- 
forces the electron-hole difference in the size of the pure 
dSC region. The range where only the dSC order pa- 
rameter is finite nearly doubles in going from U = 6t to 
U = I2t but the maximum value of Do decreases, at least 
on the hole-doped side. A stronger repulsion allows sys- 
tems with more holes to be superconducting, but at the 
same time suppresses superconductivity more effectively 
closer to half-filling. Note that the dSC order parameter 
Do should not be confused with the critical temperature: 
The maximum in Dq that appears near n = 0.9 on the 
hole-doped side does not mean that the maximum T c is 
around that doping. Instead, the maximum comes from 
the growth of Dq towards half-filling until proximity with 
the AFM phase makes it fall rapidly [22j . Finally, as seen 
in previous calculations, as U increases from 6t to 12t, the 
AF phase does not extend as far on the electron-doped 
side, [23] whereas the value of the order parameter Mq 
at half-filling increases |29j ]. 

Fig-Elshows intensity plots of the spectral functions at 
the Fermi level, for U = 8t, in the first quadrant of the 
Brillouin zone. The left illustrates a hole-doped system 
in a pure dSC phase. The spectral weight is concentrated 
along the diagonal. This is observed even without long- 
range order [26j, but is also compatible with the vanish- 
ing of the dSC gap along the diagonal. On the right, we 
display an electron-doped system in a AF+dSC phase. 
The spectral weight is depleted along the diagonal and 
concentrated near the zone boundaries ((vr, 0) and (0, 7r)). 
This is also observed in the absence of long-range order 
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FIG. 3: Intensity plot of the spectral function at the Fermi 
level, in the first quadrant of the Brillouin zone, for U — St on 
a, L — 8 cluster. Left: Hole-doped system (n = 0.93). Right: 
Electron-doped systems (n = 1.10). A Lorentzian broadening 
of 0.2t is used. 
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FIG. 4: Intensity plot of the spectral function as a function 
of uj in units of t, and wave vector, for the same parameters 
as the right-hand panel of Fig. 3. The Lorentzian broadening 
is r\ — 0.12t in the main figure and r\ = 0.044 in the inset. 



gether with the AF shadowing within the lower Hubbard 
band, are seen in Fig. 3 of Ref. [30j for n = 1.04. The 
inset, which is a higher-resolution blowup of the rectan- 
gular region shown in the figure on the (tt, it) to (tt, 0) 
segment, shows the dSC gap. 

Conclusion. This work shows, within a rigorous vari- 
ational approach that takes into account short-range dy- 
namics using clusters of different sizes, that the one- 
band Hubbard model contains the essential physics of 
the cuprates. It has dSC and AF ground states in doping 
ranges that are consistent with the observed ones for both 
electron- and hole-doping. In addition, low energy excita- 
tions manifest in the electronic spectra are shown to be 
strongly momentum dependent in essentially the same 
manner as observed using ARPES. Our results also sug- 
gest that a quantum critical point on the electron-doped 
side separates a pure dSC phase from a phase where both 
AF and dSC order parameters are non-vanishing. In the 
latter phase, we identified distinct spectral features as- 
sociated with the Mott gap and with AF and dSC long- 
range orders that should be observable by ARPES. 

We are indebted to V. Hankevych, S. Kancharla, B. 
Kyung, and M. Norman for useful discussions. This work 
was supported by NSERC (Canada), FQRNT (Quebec), 
CFI (Canada), CIAR and the Tier I Canada Research 
Chair Program (A.-M.S.T.). Computations were per- 
formed on a 200-cpu Beowulf cluster and on the 872-cpu 
Dell cluster of the RQCHP. 



|26| . but long-range order makes the plots sharper. The 
vanishing of the usual dSC gap along the diagonal is off- 
set by the AF gap. The above asymmetric behavior of 
the spectra on electron- and hole-doped sides is observed 
in Angle Resolved Photoemission Spectroscopy (ARPES) 

The Mott phenomenon, as well as the dSC and AF or- 
der parameters, each have distinct signatures in the spec- 
tral function. Fig. ^illustrates this in the AF+dSC phase 
for the electron-doped case. The Mott gap, already seen 
in ARPES 30], corresponds to an absence of states for 
all wave vectors in the interval —At < u> < —2t. The two 
main bands that disperse in the intervals —8t < uj < —At 
and — t < uj < 5t roughly correspond to those obtained 
in AF mean- field theory j^]. The AF order parameter 
manifests itself through the "shadowing" of these main 
bands by reflection about the magnetic zone boundary. 
The broad, dispersionless features around uj ~ —lQt and 
5t are remnants of the atomic limit, intimately related to 
the Mott phenomenon. The excitations occurring around 
uj ~ — 2t are absent at half- filling and from the mean-field 
solution. Their presence considerably reduces the direct 
gap at (7r/2,7r/2). Hints of these new excitations, to- 



ll] P. W. Anderson, Science 235, 1196 (1987). 

[2] E. Arrigoni et al., Phys. Rev. B 69 214519 (2004). 

[3] T. Hanaguri et al, Nature 430, 1001 (2004). 

[4] A. Paramekanti et al, Phys. Rev. B 70, 054504 (2004). 

[5] S.W. Zhang et al. 78 4486 (1997). 

[6] T.A. Maier et al, cond-mat /0409669 

[7] F. F. Assaad et al., Phys. Rev. Lett. 77, 4592 (1996). 

[8] G.-H. Gweon et al, Nature 430, 187 (2004). 

[9] M. Potthoff et al, Phys. Rev. Lett. 91, 206402 (2003). 
[10] O.Andersen et al, J. Phys. Chem. Solids 56, 1573 (1995). 
[11] D. Senechal et al, Phys. Rev. Lett. 84, 522 (2000); Phys. 

Rev. B 66, 075129 (2002). 
[12] M. Potthoff, Eur. Phys. J. B 32, 429 (2003). 
[13] A. Georges et al, Rev. Mod. Phys. 68, 13 (1996). 
[14] G. Kotliar et al, Phys. Rev. Lett. 87, 186401 (2001). 
[15] M. H. Hettler et al, Phys. Rev. B 58, R7475 (1998). 
[16] Th. Maier et aZ.,| |cond-mat /0404055 
[17] C. Dahnken et al. , |cond- mat /0309407 
[18] M. Inui et al, Phys. Rev. B 37, 2320 (1988). 
[19] E. Dernier et al, Phys. Rev. Lett. 87, 067202 (2001). 
[20] A. I. Lichtenstein et al, Phys. Rev. B 62, R9283 (2000). 
[21] G. C. Psaltakis et al, J. Phys. C 16, 3913 (1983). 
[22] B. Kyung, Phys. Rev. B 62, 9083 (2000). 
[23] Y. Dagan et al, Phys. Rev. Lett. 92, 167001 (2004). 
[24] P. Fournier et al, Phys. Rev. Lett. 81, 4720 (1998). 
[25] S.C. Zhang, Science 275, 1089 (1997). 
[26] D. Senechal et al, Phys. Rev. Lett. 92, 126401 (2004). 



■5 



[27] B. Kyung et ai, Phys. Rev. Lett. 93 147004 (2004). 
[28] J.K. Freericks et ai, Phys. Rev. Lett. 74, 186 (1995). 
[29] G. Fano et al, Phys. Rev. B 46, 1048 (1992). 
[30] N.P. Armitage et ai, Phys. Rev. Lett. 88, 257001 (2002). 



[31] F. Ronning et al, Phys. Rev. B 67, 165101 (2003). 
[32] C. Kusko et al, Phys. Rev. B 66, 140513 (2002). 



